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Abstract. Based on a simple microscopic model where the bath is in a non-equilibrium state we study 
the escape from a metastable state in the over-damped limit. Making use of Fokker-Planck-Smoluchowski 
description we derive the time dependent escape rate in the non- stationary regime in closed analytical 
form which brings on to fore a strong non-exponential kinetic of the system mode. 
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The problem of activated rate processes deals with the 
escape of a Brownian particle from a metastable state un- 
der the influence of thermal fluctuations generated by the 
immediate surrounding to which the Brownian particle is 
in close contact. Based on nonequilibrium statistical me- 
chanics Kramers [1] proposed a framework for the phe- 
nomenon which over several decades became a standard 
paradigm for theoretical and experimental investigation 
in many areas of natural science [2,3,4,5,6]. To the best 
of our knowledge majority of the post Kramers develop- 
ments of the theory have been made in the stationary 
domain with few approaches in the non-stationary regime 
within the framework of reactive flux formalism [2] . How- 
ever, few attempts have been made to deal the problem 
using nonequilibrium, non-stationary formalism [7,8,9]. 

In their work Millonas and Ray [7] proposed a the- 
oretical framework for studying the dynamics of escape 
rate from a metastable state in the over-damped limit. 
Because of the nonequilibrium fluctuations of the bath 
mode an in built fluctuating barrier appears in the effec- 
tive potential of the nonlinear Langevin equation of the 
system variable. Using path integral formalism the au- 
thors then derived a time dependent escape rate. Moti- 
vated by this work [7] we prescribe here an alternative 
method to derive the time dependent escape rate using 
the Fokker-Planck-Smoluchowski description. The object 
of the present work is twofold: First, to consider a sim- 
ple variant of system-bath model [7] to study the acti- 
vated rate processes, where the associated heat bath is in 
nonequilibrium state. The model incorporates some of the 
essential features of Langevin dynamics with a fluctuat- 
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ing barrier which has been phenomenologically proposed 
earlier [7,10]. Second, since the theories of activated rate 
processes traditionally deal with stationary bath, the non- 
stationary activated rate processes have remained largely 
overlooked so far. We specifically address this issue and 
examine the influence of initial excitation and subsequent 
relaxation of bath modes on the activation of the reac- 
tion coordinates within the framework of Fokker-Planck- 
Smoluchowski equation. In spite of the fact that our de- 
velopment bears a close kinship with the work of Millonas 
and Ray [7], it is crucial to highlight that while Millonas 
and Ray have used an explicit path integral approach to- 
wards the solution of the problem, we, on the contrary, 
implement a naive differential equation based approach 
which leads us to a closed analytical expression for the 
time- dependent escape rate. We also mention that in this 
work we have explicitly calculated the non- exponential ki- 
netics of the system mode, where the associated bath is 
not in thermal equilibrium. The closed form of the final 
expression of our approach brings with it the twin advan- 
tages of being capable of (1) handling the non-stationary 
phenomena, and (2) tracing the trajectory of how a system 
coupled with a non-equilibrium bath reaches the station- 
ary state in a computationally economic manner. 

To make the present work self consistent we describe 
the essential features of the model proposed by Millonas 
and Ray [7]. The physical situation that has been ad- 
dressed is the following. At t = 0_ , the time just before the 
system and the bath are subjected to an external excita- 
tions, the system is approximately thermalized. At i = 0, 
the excitation is switched on and the bath is thrown into 
a non-stationary state which behaves as a nonequilibrium 
bath. We follow the stochastic dynamics of the system 
mode after t > Q. The important separation of the time 
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scale of the fluctuations of the nonequilibrium bath and 
the thermal bath is that the former effectively remains sta- 
tionary on the fast correlation time scale of the thermal 
noise. 

The model consists of a system mode coupled to a 
set of relaxing modes considered as a semi-infinite dimen- 
sional system ({gfcj-subsystem) which effectively consti- 
tutes a nonequilibrium bath. This, in turn, is in contact 
with a thermally equilibrated bath. Both the baths are 
composed of two sets of harmonic oscillators characterized 
by the frequency sets {uJk} and {Oj} for the nonequilib- 
rium and equilibrium baths, respectively. The system-bath 
combination evolves under the total Hamiltonian [7] 



j k 

-x^KjQj - g{x)^qk - ^ajkQkQj- (1) 

j k j,k 

The first two terms on the right hand side describe the sys- 
tem mode. The Hamiltonian for the thermal and nonequi- 
librium baths are described by the sets {Qj, Pj} and {qj^Pj} 
for coordinate and momenta, respectively. The coupling 
terms containing kj refer to the usual system-bath linear 
coupling. The last two terms indicate the coupling of the 
nonequilibrium bath to the system and the thermal bath 
modes, respectively. In the present problem H is consid- 
ered to be classical and the temperature T is high for 
the thermally activated problem, so that the quantum ef- 
fects do not play any significant role. For simplicity we 
take mass, to = 1 in equation (1) and for the rest of the 
treatment. The form of the nonequilibrium bath, a set of 
phonons, is chosen for both simplicity and because of its 
generic relationship to many condensed matter systems. 

Elimination of equilibrium reservoir variables {Qj, Pj} 
in an appropriate way we have the equation of motion for 
the nonequilibrium bath modes as [8,11,12] 



modes with a Debye type frequency distribution of the 
form 

V{(j) = ZJ^I2ujI for \uj\<ujc 
= for \u}\ > cOc 

where lOc is the high frequency Debyc-cut-off, we finally 
arrive at the following Langevin equation of motion for 
the system mode, 

X + r{x)x + V'{x) = ^eqit) + g'{x)ineq{t). (5) 

Here r{x) is a system coordinate dependent dissipation 
constant and is given by 



r{x) =7eq -F7„eq[5'(a;)]^ 



(6) 



and S^neq (t) refers to the fluctuations of the nonequilibrium 
bath modes which effectively cause a damping of the sys- 
tem mode by an amount JneqW {x)]'^ . Equation (5) also 
includes the modification of the bare potential V{x) 



V{X) = V{X) - '^IneqQ^x). 



(7) 



equation (5) thus describes the effective dynamics of a 
particle in a modified barrier, where the metastability of 
the well originates from the dynamic coupling g{x) of the 
system mode with the noncquilibrimn bath modes. 

In order to define the dynamics described by equation 
(5) completely it is necessary to state the properties of 
the fluctuations of the nonequilibrium bath (,neq{t), which 
is assumed to be Gaussian with zero mean {£,neq(t)) = 0. 
Also the essential properties of £,neq{t) explicitly depend 
on the nonequilibrium state of the intermediate oscilla- 
tor modes {q^} through the energy density dis- 
tribution function at time t in terms of the fluctuation- 
dissipation relation for the nonequilibrium bath [7] 



qk 



+ iqk+u,lqk=9{x)+Vk{t). (2) u{.,t) = T dr{^..,,m^,,{t + r))e^' 



This takes into account the average dissipation (7) of the 
nonequilibrium bath modes qk due to their coupling to the 
thermal bath which induces fluctuations rik{t) character- 
ized by {rik{t)) = and {r]j{t)r]k{t')) = 2'ykBT5{t - t')5jk. 
In moving from equation (1) to equation (2) the cross 
terms of the form 'Y^- "fkjqj has been neglected for j ^ k. 
Proceeding similarly to eliminate the thermal bath vari- 
ables from the equation of motion of the system mode, we 
get 

X + 7egi; + V'{x) = ^eg(i) + g'ix) ^ qk, (3) 

fe 

where jeq refers to the dissipation coefficient of the system 
mode due to its coupling to the thermal bath providing 
fluctuations £,eq{t) with the properties. 



1 



ksT + e 



-7t/2 



u{u,0)--kBT 



(8) 



i^eqit)) - 0, iCeqimeqit')) = 2-f,gkBT6{t - t' 



(4) 



Now making use of the formal solution of equation (2) 
which takes into account of the relaxation of the nonequi- 
librium modes, and integrating over the nonequilibrium 



where 0) — ^fc^r] is a measure of departure of en- 

ergy density from thermal average aX t — Q. The expo- 
nential term implies that deviation due to the initial ex- 
citation decays asymptotically to zero as t ^ 00, so that 
one recovers the usual fluctuation-dissipation relation for 
the thermal bath. With the above specification of corre- 
lation function of ^neq, equation (8) thus attributes the 
non-stationary character of {gfc}-subsystem. 

On time scales larger than the inverse friction coef- 
ficient l/7e(j, we can in most particular cases consider 
the over-damped limit of the Langevin equation. This in 
turn corresponds to the adiabatic elimination of the fast 
variables, inertia term, from the equation of motion by 
putting X = for homogeneous systems. In contrast, for 
the case of inhomogeneous system the above method of 
elimination does not work properly and Sancho et al [13] 
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have given a proper prescription for the ehmination of fast 
variables. Using the method of Sancho et al the formal 
master equation for the probability density of the process 
P{x,t) = {p{x,t)) is given by 



dP 

'dt 




1 d 



r{x) dx r{x) 

g'jx) d g'{x) 

r{x) dx r{x) 

9'{x)g"{x) 



P 



r^x) 



(9) 



where r = {[u{ui 0, 0)/2kBT] — 1} is a measure of the 
deviation from equilibrium at the initial instant, equation 
(9) is the Fokker-Planck-Smoluchowski equation where the 
associated bath is in nonequilibrium state, and is the first 

key result of this paper. Under stationary condition (at 
t — > oo) dP/dt = and the stationary distribution obeys 
the equation 



kBT^^^ + V\x)P,,{x)=Q 



which has the solution 

P,t(x) =7Vexp 



keT 



V'{x')dx' 



(10) 



(11) 



where N is the normalization constant. In ordinary Strato- 
novich description the Langevin equation corresponding 
to the Fokker-Planck-Smoluchowski equation (9) is given 

by 



X = — 



V'{x) D{t)g"{x)g'{x) 



r{x) 



r^{x) 



+ 



1 



r{x) 



r{x) 



where 



(12) 

(13) 



is the time-dependent diffusion constant due to the relax- 
ation of nonequilibrium bath. Using equations (9) or (12) 
the escape rate from a metastable state can be calculated 
via steepest descent method [14] 



(14) 



where k is the Kramers activation rate, with Ei, = V'{xi,) 
V'{xo) is the modified activation energy and (Df, = [V""(a;b)]^ 
ijjo = [V"{xo)]^^^ are the modified frequencies at the bar- 
rier top and the bottom of the potential well, respectively. 
Xb denotes the position of the barrier top and xq is the 
position of the bottom of the potential well. In equation 
(14) r has been evaluated at the top of the barrier. In 
the absence of the nonequilibrium bath (14) reduces to 
standard Kramers' expression [1], 



To obtain the time dependent rate k{t), let us consider 
that the time dependent solution of equation (9) is given 
by 

P{x,t) = Pst{x)e-*^'^ (16) 

where (/) is a function oft only and limi^oo (t'i't) = 0. Psi{x) 
is the steady state solution of equation (9). Substitution 
of (16) in (9) separates the space and time parts and we 
have the equation for (j){t) as 



— ^^'^^^^ — constant = a{say) 
which after integration over time gives 



^e-7t/2 



(17) 



(18) 



where a can be determined by initial condition. The time 
dependent solution of equation (9) therefore reads 



P{x,t) = Pst{x)exp 



2a 
( 

7 



-7f/2 



(19) 



To determine a we now assume that [8] just at the mo- 
ment the system (and the non-thermal bath) is subjected 
to external excitation at t = 0, the distribution must co- 
incide with the usual Boltzmann distribution where the 
energy term in the Boltzmann factor in addition to the 
usual kinetic and potential energy terms, contains the ini- 
tial fluctuations of energy density Au[= u{u),Q) — ^ksT] 
due to excitation of the system at t = 0. This gives a = 
("f /2) (Au/ kgT) . a is thus determined in terms of relaxing 
mode parameters and fluctuations of the energy density 
distribution at t = 0. The time dependent rate is then 
derived as [14] 



k{t) = k exp 



Au 
k^' 



(20) 



27r7 



exp 



Eb \ 
'ksT) ■ 



(15) 



where Au is the measure of the initial departure from the 
average energy density distribution due to the preparation 
of the non-stationary state of the intermediate bath modes 
as a result of excitation at i = and k is given by equation 
(14). The above result, equation (20), which is the second 
key result of this paper, illustrates a strong nonexponential 
kinetic of the system mode undergoing a non-stationary 
activated rate processes in the over-damped regime. The 
origin of this is an initial preparation of nonequilibrium 
mode density distribution which eventually relaxes to an 
equilibrium distribution, equation (20) implies that the 
/^i,nitial transient rate is different from the asymptotic steady 
state Kramers' rate. The sign of Au determines whether 
the initial rate will be faster or slower than the steady state 
rate. This is because there exists a time lag for the non- 
thermal energy gained by the few nonequilibrium modes 
by sudden excitation to be distributed over a range before 
it becomes available to the reaction coordinate as thermal 
energy for activation. 

In conclusion, based on a system-reservoir model, where 
the reservoir is in a non-equilibrium state, we have pro- 
vided an analytic model to derive the closed time-dependent 
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escape rate from a metastable state induced by non-equilib- 
rium fluctuations. We have explicitly calculated the non- 
exponential kinetics of the system mode, where the associ- 
ated bath is not in thermal equilibrium. Our methodology 
takes care of the non-stationary phenomena, and simul- 
taneously traces the barrier dynamics of a system when 
it is coupled with a non-equilibrium bath. Not only that 
our approach may serve as a potential avenue towards 
the explanation of non-stationary transport processes and 
rachet problems envisaged in various chemically and bio- 
logically interesting systems. The work in this direction is 
in progress in our group. 
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